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We report results for the interaction measure, pressure and energy density for nonzero temperature 
QCD with 2+1 flavors of improved staggered quarks. In our simulations we use a Symanzik 
improved gauge action and the Asqtad O (a 2 ) improved staggered quark action for lattices with 
temporal extent N t = 4 and 6. The heavy quark mass m s is fixed at approximately the physical 
strange quark mass and the two degenerate light quarks have masses m u( i = 0.1m s or 0.2m s . The 
calculation of the thermodynamic observables employs the integral method where energy density 
and pressure are obtained by integration over the interaction measure. 
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1. Introduction 

The equation of state (EOS) is important for phenomenological models of quark-gluon plasma 
formation and decay, which is currently under experimental study at RHIC and elsewhere. We 
have determined the EOS with the Asqtad quark action [jl|] for 2+1 flavors, combined with a one- 
loop Symanzik improved gauge action [gj. The Asqtad action is well suited for high temperature 
studies since it has excellent scaling properties and much better dispersion relations in the free case 
than the standard Wilson or staggered actions, which means decreased lattice artifacts above the 
transition. 

Our nonzero temperature studies are at N t = 4 and 6. Even with Asqtad improvement, for 
T < T c , an ,/V, = 4 lattice has a badly split pion taste multiplet with some members heavier than the 
kaon, which makes for questionable strange quark physics. At N t = 6 the taste-splitting is about 
half as large. One of our goals was to determine to what extent the increase in N t from 4 to 6 
influences the EOS. 

2. Action 

The fermion part of the action we use is effectively written as: 

S f = -Y,(nf/4)Ti\n[M(am f ,U,u )], (2.1) 

/ 

where M(amf,U,uo) is the fermion matrix corresponding to the Asqtad2+1 flavor staggered action. 
The gauge part is defined as: 

Sg = P E (l-^v)+At E (l-fl M v)+frh E (l-C M va). (2-2) 

x,ji<v x,fi<v x,ii<v<a 

The gauge couplings above are p = 10/g 2 , fa = -0(1 +0.4805oO/(20w§), j8 ch = -0.03325)3 a s /u 2 Q , 
with a s = -41n(>o)/3.0684 and u = (P) l/A . 

For our simulations we use the dynamical R-algorithm [||] with step-size equal to the smaller of 
0.02 and 2am uc i/3. Our aim is to generate zero and nonzero temperature ensembles of lattices with 
action parameters chosen so that a constant physics trajectory {m n jm p = const) is approximated. 
Along the trajectory the heavy quark mass is fixed close to the strange quark mass. We work with 
two such trajectories: m uc j = 0.2m s (m n /m p 0.4) and m MC / = 0.\m s {m^/rrip « 0.3). 

3. Parameterization of the Constant Physics Trajectories and Run Parameters 

Our trajectories are intended to approximate constant zero temperature physics. The con- 
struction of each trajectory begins with "anchor points" in j8, where the hadron spectrum has been 
previously studied [Q] and the lattice strange quark mass has been tuned to approximate the correct 
strange hadron spectrum. We adjusted the value of am uc i at the anchor points to give a constant 
(unphysical) ratio m K /m p . Between these points the trajectory is then interpolated, using a one- 
loop renormalization group inspired formula. That is, we interpolate \n(am s ) and \n{am U( j) linearly 
in p. Since we have three anchor points for the m uc j = 0.2m s trajectory, namely /3 = 6.467, 6.76, 
and 7.092, our interpolation is piecewise linear. For the trajectory m ui ] = 0.\m s we use two anchor 
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/3 


am ud 


am s 


«o 


i -f-\> 


Vr=o 


*6.300 


0.0225 


0.1089 


0.8455 


12 3 x 6 


12 4 


*6.350 


0.0206 


0.1001 


0.8486 


12 3 x 6 


12 4 


6.400 


0.01886 


0.0919 


0.8512 


12 3 x 6 




6.433 


0.01780 


0.0870 


0.8530 


12 3 x 6 




*6.467 


01676 


0.0821 


0.8549 


16 3 x 6 


16 3 x 48 


6.500 


01580 


0.0776 


0.8568 


12 3 x 6 




*6.525 


01510 


0.0744 


8580 


12 3 x 6 


12 4 


6.550 


0.01450 


0.0713 


0.8592 


12 3 x 6 




*6.575 


0.01390 


0.0684 


0.8603 


12 3 x6 


16 4 


6.600 


0.01330 


0.0655 


0.8614 


12 3 x6 




*6.650 


0.01210 


0.0602 


0.8634 


12 3 x6 


20 4 


6.700 


0.01110 


0.0553 


0.8655 


12 3 x6 




*6.760 


0.01000 


0.0500 


0.8677 


20 3 x6 


20 3 x 64 


7.092 


0.00673 


0.0310 


0.8781 


12 3 x6 




7.090 


0.00620 


0.0310 


0.8782 




28 3 x 96 



Table 1: Run parameters of the trajectory with »y = 0.2m s at N t = 6. The asterisk indicates parameter sets 
for which both zero and nonzero temperature runs were performed. 



p 


am ud 


am s 


«0 


V T^0 


VV = o 


*6.300 


0.01090 


0.1092 


0.8459 


12 3 x6 


12 4 


*6.350 


0.00996 


0.0996 


0.8491 


12 3 x6 


12 4 


6.400 


0.00909 


0.0909 


0.8520 


12 3 x6 




*6.458 


0.00820 


0.0820 


0.8549 


16 3 x6 


12 4 


6.500 


0.00765 


0.0765 


0.8570 


12 3 x6 




*6.550 


0.00705 


0.0705 


0.8593 


12 3 x6 


20 4 


6.600 


0.00650 


0.0650 


0.8616 


12 3 x6 




*6.650 


0.00599 


0.0599 


0.8636 


12 3 x6 


24 4 


6.700 


0.00552 


0.0552 


0.8657 


12 3 x6 




*6.760 


0.00500 


0.0500 


0.8678 


20 3 x6 


24 3 x 64 


*7.080 


0.00310 


0.0310 


0.8779 


18 3 x6 


40 3 x 96 



Table 2: Run parameters of the trajectory with m uc \ = 0.lm s at N t = 6. The asterisk indicates parameter sets 
for which both zero and nonzero temperature runs were performed. 



points at jS = 6.458 and 6.76. For both trajectories, for values of j3 out of the interpolation intervals, 
the parameterization formulas are used as extrapolations appropriately. The run parameters of the 
two trajectories at different N t are summarized in Tables 1, 2 and 3. 



4. Integral Method for the EOS derivation 

To derive the analytic form of the EOS we employ the integral method [§]. We start from the 
thermodynamic identities: 



£V = - 



dlnZ 



d(l/T) 



p dlnZ 
T 



dV 



InZ 



I = e-3p=- 



T dlnZ 



V <?lna' 



(4.1) 
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p 


am uc j 


am s 


«o 


i -f-\> 


Vr=o 


*6.000 


01980 


0.1976 


0.8250 


12 3 x 4 


12 4 


*6.050 


0.01780 


0.1783 


0.8282 


12 3 x 4 


12 4 


6.075 


0.01690 


0.1695 


0.8301 


12 3 x 4 




*6.100 


0.01610 


0.1611 


0.8320 


12 3 x 4 


12 4 


6.125 


01530 


0.1533 


8338 


12 3 x 4 




*6.150 


01460 

V.' . \J 1 IV' V.' 


0.1458 


8356 


12 3 x 4 


12 4 


6.175 


01390 


0.1388 


0.8374 


12 3 x 4 




*6.200 


01320 


0.1322 


0.8391 


12 3 x 4 


12 4 


6.225 


01260 


0.1260 


0.8407 


12 3 x 4 




*6.250 


0.01200 


0.1201 


0.8424 


12 3 x 4 


12 4 


6.275 


0.01140 


0.1145 


0.8442 


12 3 x 4 




*6.300 


01090 


0.1092 


0.8459 


12 3 x 4 


12 4 


*6.350 


00996 

v ' . \ / v./ y y v/ 


0.0996 


0.8491 


12 3 x 4 


12 4 


6.400 


0.00909 


0.0909 


0.8520 


12 3 x 4 




*6.458 


0.00820 


0.0820 


0.8549 


12 3 x4 


12 4 


6.500 


0.00765 


0.0765 


0.8570 


12 3 x4 




*6.550 


0.00705 


0.0705 


0.8593 


12 3 x4 


20 4 


6.600 


0.00650 


0.0650 


0.8616 


12 3 x4 




*6.650 


0.00599 


0.0599 


0.8636 


12 3 x4 


24 4 


6.700 


0.00552 


0.0552 


0.8657 


12 3 x4 




*6.760 


0.00500 


0.0500 


0.8678 


12 3 x4 


24 3 x 64 



Table 3: Run parameters of the trajectory with m ut j = 0.1m, at N t = 4. The asterisk indicates parameter sets 
for which both zero and nonzero temperature runs were performed. 

where the derivative with respect to In a is taken at constant m n /m p and the partition function is Z. 
Using the above identities and the explicit form of Z we obtain: 



la' 



. 6 ^_A{P)-12fLA{J?}-A{C) 



dlna 



pa 



n f 



lna 



d{mfa) 
dlna 

5 dp 
dlna 



dlna 



dlna 

duo, I _ dM 

A( v/ y/ 

dlna \ duo 



(4.2) 



A (P) + 12^A (R) + 16^A (C) 



dlna 



rtf 



f 



d(nifa) 



dlna 



A (WW) 



dlna 

duo / _ dM 
^ dlna \ duo 



dlna', 



ea 4 = (I+3p)a 4 



(4.3) 



(4.4) 



In the above expressions the symbol A in front of the lattice observables stands for the difference 
between the values of those observables at nonzero and zero temperatures. In the pressure expres- 
sion, the lower integration endpoint lnao is set where the zero temperature subtracted value of la 4 
is zero within errors at coarse lattice spacings. To calculate the EOS, in addition to the lattice glu- 
onic and fermionic observables in the above analytic forms, we need to calculate the derivatives 
dp p ]/dlna, dfi It /dlna, dp c h/d\na, dinifd) /dlna and du /dlna. For this purpose we take deriva- 
tives of the ln(ara H£ /) and ln(ara v ) trajectory parameterizations, polynomial fits to «o(j3) for both 
trajectories, and the updated version of the ln(ri/a) fitting formula in [Q] shown below: 



ln(n/a) = Coo + Cio(/3 -7)+C 2 i(2am„ rf + flm. v )+C 2 o(j3 -7) 2 



(4.5) 
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where n = 0.317(7) (3) fm, C 00 = 1.261(3), C w = 0.939(8), C 2 \ = -0.86(3) and C 20 = -0.25(2). 
The above fit (4.5) giving the relative lattice scale is based on measurements of the static quark 
potential at zero temperature for a large set of /3 and quark masses. The absolute scale is fixed from 
a determination of the bottomonium spectrum The fit has % 2 /DOF 1.1. 

5. EOS results 

Figure 1 summarizes our results for the EOS. The errors on all data points are calculated using 
the jackknife method, and we ignore insignificant errors on the derivatives of the bare parameters 
with respect to the lattice scale discussed at the end of the previous section. For the points where 
there is no zero temperature run, local interpolations are made to calculate the zero temperature 
corrections to the interaction measure. The integration of the interaction measure to obtain the 
pressure is done using the trapezoid method. 

The comparison between the N t = 4 and 6 cases for the m u d = 0. bn s trajectory shows that there 
is not a significant difference between them, except in the interaction measure near the transition 
region. The EOS results from the two different physics trajectories are very similar and for the 
temperature interval we studied, the deviation from the 3 flavor Stefan-Boltzmann values are large. 

In the temperature region where we have data, we consider our pressure results to be in general 
agreement to a previous p4-action N t = 4 calculation [Qj. 

6. Conclusions 

We have calculated the EOS for 2+ 1 dynamical flavors of improved staggered quarks (m uc i / m s = 
0.1 and 0.2) along trajectories of constant physics, at N t = 4 and 6, where the latter is the first result 
of its kind. Our results show that the N t = 4 and N t = 6 results are quite similar except in the 
crossover region where the interaction measure is a bit higher on the finer N t = 6 lattice. We also 
do not see significant differences between the EOS results from the two physics trajectories. We 
find large deviations from the 3 flavor Stefan-Boltzmann limit. Our results are comparable with 
previous calculations [f7|]. 
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Figure 1: The EOS with 2+1 dynamical quark flavors for both constant physics trajectories. The energy 
density, pressure and the interaction measure symbols are diamonds, squares and circles, respectively. The 
data at N t = 6 is in black and at = 4 in red. All data with filled symbols belong to the m uc i = 0.2m, 
trajectory and the empty symbols belong to the m U£ / = O.lrn, trajectory. The Stefan-Boltzmann lattice limits 
for 3 flavors of quarks at N t =4 and 6 are represented by solid lines - red and black, respectively. The black 
dashed line is the continuum Stefan-Boltzmann limit for 3 flavors of quarks. 
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